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ABSTRACT 

The three-point correlation function (3PCF) can now be measured in large galaxy 
redshift surveys, but in three dimensions its interpretation is complicated by the pres- 
ence of redshift-space distortions. I investigate the projected 3PCF, where these distor- 
tions are eliminated by integrating over the redshift dimension, as is commonly done 
for the two-point correlation function. The calculation of the projected 3PCF from the 
real-space, three-dimensional bispectrum is greatly simplified by expanding both quan- 
tities in Fourier components, analogous to Szapudi's expansion of the three-dimensional 
quantities in multipole components. In the weakly nonlinear regime, the bispectrum 
can be well represented by the first few Fourier components. There is a well-known rela- 
tion between the reduced 3PCFs of matter and galaxies in the weakly nonlinear regime, 
which can be used to infer galaxy bias factors if the real-space three-dimensional galaxy 
correlation functions (two-point and three-point) can be measured. I show that the same 
relation holds for the reduced projected 3PCFs if these are properly defined. These re- 
sults should aid determinations of galaxy bias from large redshift surveys by eliminating 
the complication of redshift-space distortions. 

Subject headings: cosmology: theory — dark matter — galaxies: formation — clustering 
— galaxies: statistics — galaxies: halos - large-scale structure of universe 



1. Introduction 

The three-point correlation function (3PCF), or its Fourier transform, the bispectrum, is a 
valuable complement to two-point statistics in characterizing galaxy clustering. The behavior of 
the 3PCF of the matter is well understood in perturbation theory (e.g., Fry 1984; Goroff et al. 1986; 
Bernardeau 1992; Jain & Bertschinger 1994; Scoccimarro & Frieman 1996; Scoccimarro et al. 1998; 
Scoccimarro &: Frieman 1999; Bernardeau et al. 2002). For Gaussian initial conditions, second- 
order perturbation theory predicts that the amplitude of the 3PCF scales like the square of the 
amplitude of the two-point correlation function (2PCF) £ (Peebles 1980), and this scaling is one of 
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the fundamental tests of the Gaussian primordial fluctuations predicted by inflationary cosmology. 
Nonlinear gravitational evolution produces anisotropic, filamentary structures, so elongated triangle 
configurations have higher amplitude in the three-point statistics. Galaxy bias, a difference between 
the distributions of galaxies and matter, can alter three-point statistics, but it tends to boost or 
suppress amplitudes for all triangle configurations equally, at least on large scales. The triangle 
shape dependence of the 3PCF or bispectrum, in combination with the 2PCF or power spectrum, 
therefore becomes a diagnostic for galaxy bias and the matter clustering amplitude (Fry 1994). 

Early measurements of the 3PCF or bispectrum were based on angular catalogs (e.g., Peebles & 
Groth 1975; Fry & Peebles 1980; Jing, Borner, & Mo 1991; Fry 1994). Galaxy redshift surveys like 
the Two-Degree Field Galaxy Redshift Survey (2dFGRS; Colless et al. 2001) and the Sloan Digital 
Sky Survey (SDSS; York et al. 2000), are now large enough to allow measurements of the redshift- 
space 3PCF or bispectrum with high signal-to-noise ratio (e.g., Verde et al. 2002; Jing & Borner 
2004; Kayo et al. 2004). However, this brings in the additional complication of distortion by the 
peculiar motions of galaxies. Even for the 2PCF or power spectrum, the effect of nonlinear redshift- 
space distortions persist to remarkably large scales (Cole, Fisher, & Weinberg 1994). Distortions 
of the 3PCF or bispectrum are more complex, and while some models of these distortions exist 
(e.g., Scoccimarro et al. 2001; Verde et al. 2002), it is not clear that they are accurate at the level 
of precision afforded by current data. An obvious way to circumvent redshift-space distortions 
is to project the 3PCF over the redshift direction, as is commonly done for the projected 2PCF 
w p (r p ) (Davis & Peebles 1983). We note that projection of redshift surveys is not the same as 
simply measuring from the parent angular catalog — the use of galaxy redshifts to obtain physical 
projected separations greatly reduces the noise and yields a quantity more closely related to the 
three-dimensional 3PCF. Projected 3PCFs have been measured for the Las Campanas Redshift 
Survey (LCRS) by Jing & B6rner (1998) and for the 2dFGRS (Jing k Borner 2004) on relatively 
small scales. By comparing the measurements with the predicted 3PCFs of matter in iV-body 
simulations, they find that the observed galaxy 3PCFs are lower than the predicted matter 3PCFs, 
indicating the need of galaxy bias to explain the data. 

Since the three-dimensional 3PCF is already a complicated object, projection may seem to forgo 
any chance of analytic treatment. However, I show here that the relation of the projected 3PCF to 
the (three-dimensional) bispectrum is relatively straightforward if the bispectrum is decomposed 
into Fourier moments (§ 2). My analysis is directly analogous to that of Szapudi (2004), who 
introduces multipole expansion of the three-dimensional 3PCF and bispectrum and shows it to 
be useful in simplifying the relation between them, in characterizing the triangle configuration 
dependence, and in constraining galaxy bias. The spirit of this paper parallels that of Szapudi 
(2004). In § 3 I show that bias effects remain simple for projected 2PCFs and 3PCFs on large 
scales, which makes it attractive to constrain galaxy bias using projected quantities. I briefly 
summarize my results in § 4. 
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2. Fourier Expansion and Projected 3PCFs 



Before going into details, I emphasize that the meaning of projection in this paper has subtle 
differences from the usual angular correlation. Our goal is to describe the projection of the measured 
redshift-space correlation functions from a galaxy redshift survey and to use it to study problems 
like galaxy bias. Since galaxies with different properties cluster differently (e.g., red galaxies are in 
general more strongly clustered than blue galaxies) , it is desirable to construct volume- limited sam- 
ples that uniformly represent galaxies of a given type. By doing so, we can infer bias information as 
a function of galaxy type, rather than some average over all types of galaxies (weighted in a com- 
plicated way) like that from flux-limited galaxy samples. This approach helps the understanding 
of galaxy bias and aids comparisons with theoretical structure formation models. Redshift-space 
correlation functions can be measured by comparing the galaxy distribution with a distribution of 
random points that has the same geometry and selection function as the galaxies (e.g., Landy & 
Szalay 1993). The projected correlation function can be formed by integrating the redshift-space 
correlation function along the redshift direction. The selection function does not enter the inte- 
gration because its effect is included in the error budget of the measured redshift-space correlation 
function. In contrast, one needs to take into account the selection function in modeling angular 
correlation measurements of galaxy clustering. Nevertheless, our results in this paper can be easily 
generalized to include selection function to model angular clustering (see the discussion in § 3 and 
Fry & Thomas 1999). 

Since redshift-space distortions conserve numbers of pairs and triplets, 2PCFs and 3PCFs pro- 
jected from real-space three-dimensional correlation functions are identical to those from redshift- 
space three-dimensional correlation functions. This vastly simplifies theoretical predictions of pro- 
jected 2PCFs and 3PCFs — we only need to project the real-space three-dimensional correlation 
functions, without calculating redshift-space distortions. In practice, the radial extent of the pro- 
jection may not be infinite, but it should be large enough to minimize any residual redshift-space 
distortions. We generally assume ideal infinite projections and plane-parallel geometry (the distant 
observer approximation) in our derivation. 



We start by reviewing the calculation of the projected 2PCF and its relation to the fluctuation 
power spectrum. 

The (three-dimensional) 2PCF £(r) is the Fourier transform of the fluctuation power spectrum 



where r is the pair separation. Equation (1) holds for real-space quantities, which are isotropic, 
or for redshift-space quantities, which are not. Both r and k can be decomposed into components 



2.1. Projected 2PCFs 



P(k) 




(1) 
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perpendicular and parallel to the line of sight: r = r p + ry and k = k p + ky. The projected 2PCF 
is obtained by integrating £(r) along the line of sight, 

/+oo 
dr||C(r p + r||). (2) 

After substituting equation (1) into equation (2) and putting variables into the form of perpendic- 
ular and parallel components, we have 

i r r+ca r+oo 

w p (r p ) = —JdX^ rp J dfe N ^(kp + kn) y dr\\ e l Vll (3) 

= jhfj d2k p p t k P y krrp - ( 4 ) 

Equation (4) follows from equation (3) because the rightmost integral in equation (3) is just 
2tt times the Dirac (5-function 5o(k^). Equation (4) states that the projected 2PCF is the two- 
dimensional Fourier transform of the power spectrum, and because we have projected out redshift- 
space distortions, we can use the isotropic, real-space P(k) in the integral. We can evaluate 
equation (4) in polar coordinates. The angular part can be calculated by using the expansion 
exp(ik-r) = Y2n=-<x> J n (kr)i n exp[in(<f> — <3?)] (plane waves in terms of cylindrical waves), where 
(j) — <3? is the angle between k and r and J n (x) is the Bessel function of integer order. Finally, 
equation (4) reduces to a one-dimensional integral, 

w p (r p ) = J —dkP(k)J (kr p ). (5) 

This equation mimics the relation for the real-space, three-dimensional 2PCF which involves the 
spherical Bessel function jo(x). This kind of result can be found in papers that deal with projected 
observations, such as variants of Limber's equation (e.g., Baugh & Efstathiou 1993). 



2.2. Projected 3PCFs 

The projected 3PCFs can be derived in a similar way to the 2PCFs. 

The three-dimensional 3PCF is the Fourier transform of the bispectrum -B(ki, k2, k3), 

e [3l (ri,r 2 ,r 3 ) = JL J d^Ad 3 ^ B(k u k 2 , k 3 )e^ kl - ri+k2 - r2+k3 - r3 )^(k 1 + k 2 + k 3 ). (6) 

The Dirac 5-function selects wavevectors (ki,k2,k 3 ) that form a triangle. From now on, we adopt 
the notation that a subscript of one number represents one side for a wavevector triangle, while 
it represents one vertex for a real-space triangle. Each side of a real-space triangle is denoted by 
subscript of two numbers, e.g., = i-j — rj = 1, 2, 3). Equation (6) reduces to 

e [3l (ri,r 2 ,r 3 ) = ^ J d^^B^^y^^^, (7) 
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where k3 does not appear explicitly in the bispectrum because the wavevector triangle is fully 
determined by ki and k2 (ki + k2 + k3 = 0). 

Similar to the projected 2PCF, the projected 3PCF for a projected triangle with (r pl , r p2 , r p3 ) 
as vertices can be obtained by integrating ^ along the line of sight. One vertex (e.g., r p3 ) can be 
fixed. We then have the following expression for the projected 3PCF, 

/+oo f+OO 
dr\\! dr|| 2 ^ 3 ](r p i +r|| 1 ,r p2 +r,| 2 ,rp 3 ). (8) 

-oo J — oo 

Substituting equation (7) into the above, changing to cylindrical coordinates, and following the 
procedures used to go from equation (3) to equation (4), we can derive the relation between the 
projected 3PCF and the bispectrum 

(r pl , r p2 , r p3 ) = J d 2 k pl d 2 k p2 B(k pl , k p2 )e ik ^ W k ^ 23 . (9) 

That is, the projected 3PCF is the two-dimensional Fourier transform of the bispectrum. Equa- 
tion (9) is a four-dimensional integral, and it is clearly desirable to simplify the calculation, by 
reducing its dimensionality. Note that B(k p i,k p2 ) is just the usual three-dimensional, real-space 
bispectrum evaluated at ki = k p i, k2 = k P 2, and k 3 = — (k p i + k P 2). 

Szapudi (2004) introduces multipole expansion of the three-dimensional three-point statistics 
and shows that the 3PCF can then be put in a simple form (see also Verde, Heavens, & Matarrese 
2000 for a similar expansion in the bispectrum). Consequently, the whole calculation becomes 
simple because only a few multipoles are needed for accurate convergence. In the same spirit 
as Szapudi (2004), we introduce Fourier expansion of the projected three-point statistics. The 
bispectrum can be expanded as 

+oo 

B(k 1 ,k 2 , ( f>)= ^ B n (k u k 2 )e in t (10) 

n=— oo 

where <j> is the angle between ki and k2, and the coefficient B n (ki, k 2 ) can be obtained through 

B n (h,k 2 ) = -!- 1^ dcpB{k 1 MAV n<t> = ^ 1^ d<f>B(k 1 ,k 2 ,(f>) cos n<p. (11) 
2vr J 2tt J 

We have used the fact that the bispectrum has the symmetry B(k±, k 2 , — 4>) = B(ki,k 2 ,cf)), so 
B_ n = B n and the expansion is a cosine Fourier expansion. 

We rewrite equation (9) in polar variables, replace B with its Fourier expansion, and expand 
the two exponentials in the same way as we do in deriving equation (5) (it is useful to write the 
angle between two vectors as the difference of their polar angles). After integrating the angular 
part using the orthogonality of exp(in0), we find that 

+oo 

wW(r P i3,r p23 ,<£)= 4 3] >pi3^2 3 y"*, (12) 

n=— oo 
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where $> is the angle between r p i3 and r p 23, and the Fourier coefficient 

roc roc ^ 

Wpjn( r pi3,r p 23) = J ^ dkl J ^ dk 2(-l) n B n (ki,k 2 )Jn{kir pl3 )J n (k 2 r p 23). (13) 

Equation (13) is in a form that resembles the projected 2PCF in equation (5). The essence here 
is that, by Fourier expansion of the bispectrum B and cylindrical wave expansion of the plane 
wave exp(ik p -r p ), the angular dependence is separated and integrated. There are also similarities 
between the above equation and equation (3) in Szapudi (2004) for the multipole coefficient of the 
three-dimensional 3PCF. We see that the calculation of the projected 3PCF from a given three- 
dimensional bispectrum requires essentially the same amount of work as the three-dimensional 
3PCF (and a similar remark holds for 2PCF and the three-dimensional power spectrum). For 
3PCFs, what we need to calculate are expansion coefficients, where only one-dimensional (for 
expansion coefficients of the bispectrum) and two-dimensional (for coefficients of the correlation 
function) integrals are involved. I note that the same result (equation [13]) is also obtained inde- 
pendently by I. Szapudi (2004, private communication). 

If the effect of finite projection has to be taken into account, the Dirac 5-function we use 
to derive equations (4) and (9) is replaced by r\\ max jo(k\\r\\ max )/ir, if the projection is performed 
for line-of-sight separations from —r^ max to ry >max . Equations (4) and (9) become three- and 
six-dimensional integrals, respectively. The Fourier expansion can still be adopted to reduce the 
calculation of the projected 3PCF to four-dimensional integrals. 



As an example of the advantage of the Fourier expansion, Figure 1 shows the first few Fourier 
expansion coefficients of the matter bispectrum for k\ = 2k 2 = 0.05/iMpc -1 and the bispectrum 
recovered from using Fourier components up to n = 3, 5, and 10 in the weakly nonlinear regime. 
The matter bispectrum in this regime can be written in the following form (e.g., Fry 1984), 



B(k 1 ,k 2 ,(f>) 



3 1 
2~2^ 



+ 



h. + -i 

k 2 k\ 



cos <f> + -(1 — ju) cos 2(p 



P(k 1 )P(k 2 )+ perm., (14) 



where P(k) is the linear matter power spectrum and fj, reflects the weak dependence on cosmology 
(H = 3tt m 1/M0 /7 » 3/7 for a spatially flat universe; Kamionkowski & Buchalter 1999, see also 
Matsubara 1995). The first term has Fourier components only up to n = 2. The other permutations 
have higher frequency components because k$, ki-k 3 , and k2-k3 are all functions of cos</>. 



The results in Figure 1 are analogous to those of Szapudi (2004), where the bispectrum is 
viewed as contributions by multipoles. In fact, the nth-order Legendre polynomial has Fourier 
components up to n. In Figure 1 we have adopted a linear power spectrum in the parameterization 
of Efstathiou, Bond, & White (1992) with a primordial fluctuation power-law index n s = 1 and 
a shape parameter V = 0.21. We see that only a few, low-frequency Fourier components are 
significant. For example, truncating the Fourier series up to n = 5 can recover the bispectrum with 
a fractional error of ~1% for k\ = 2k 2 = 0.05/iMpc _1 . In the weakly nonlinear regime, the galaxy 
bispectrum can be obtained by a linear combination of the matter bispectrum and the quantity 
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P(ki)P(k2) + P(k2)P(ks) + P(ks)P(ki) with combination coefficients depending on galaxy bias 
(e.g., Fry 1994). In modeling galaxy 3PCFs, it is useful to expand this quantity in Fourier series, 
too. Again, the first few Fourier coefficients are enough. 



3. Reduced Projected 3PCFs and Galaxy Bias 

The reduced 3PCF in real space is defined as 

£[ 3 l(ri,r 2 ,r 3 ) 

Q(ri ' r2 ' r3) = e(ri 2 )e(r 2 3)+^(r 2 3)e(r 3 i) + e(r3i)e(ri 2 )- (15) 

For galaxies, this quantity has contributions from both gravity and galaxy bias (e.g., Fry & 
Gaztahaga 1994; Fry 1994; Juszkiewicz et al. 1995). For Gaussian initial conditions, in the weakly 
nonlinear regime, if the local galaxy overdensity 5 g is expanded in terms of the local matter over- 
density 5 m as S g = J2bn5m/ n * (local bias model), the reduced 3PCFs of galaxies and matter have 
the relation 

1 b 2 
Q g (ri,r 2 ,r 3 ) = -Q m (ri, r 2 , r 3 ) + ^, (16) 

where b = b\ is the linear galaxy bias factor and 6 2 is the lowest-order nonlinear bias factor. In 
redshift space, the 3PCF of galaxies, which is easy to measure, differs from the real-space one 
because of the additional contribution from peculiar motions of galaxies, and the above relation 
no longer holds. Constraints on galaxy bias from observations in redshift space therefore require 
accurate modeling of the peculiar motion of galaxies. Since projection along the line of sight 
essentially eliminates the effect of redshift-space distortions, it is interesting to see whether we can 
infer b and 6 2 using measurements of projected quantities. 

At first glance, it seems impossible to directly perform projections on both sides of equa- 
tion (16), since the 6 2 /6 2 term leads to divergence. Furthermore, the projection of Q from redshift 
space is not, in general, the same as from real space. We therefore project the 2PCFs and 3PCFs 
themselves, before taking ratios. To do so, we substitute the definition of Q (equation [15]) into 
equation (16) and multiply both sides by [^ 9 (r p i 2 )^ ff (r p23 ) + perm.]. On large scales, by using the 
relation £ 9 = 6 2 £ m , we find 

4 3 ]( ri ,r 2 ,r 3 ) = 6 3 ^( ri ,r 2 ,r 3 ) + | &(n 2 )£ fl (r 23 ) + ^(r 23 )^(r 31 ) + ^(r 31 )^(r 12 )] . (17) 

Now we can perform projection on both sides by fixing r p3 , as we do in equation (8). After 
the projection, £l 3 l becomes toj, 3 ' and it is easy to show that the product Cg{ r ij)Cg{ r jk) becomes 
Wp,g{rpij)wp^g{r p jk). We then divide both sides by [u> PiS (r p i 2 )u;p j9 (rp 23 ) + perm.] and use the linear 
bias relation w p>g = b 2 w p>m on large scales. The result is just equation (16) in terms of the reduced 
projected 3PCF Q p , where Q p is formed as 

g( r u w§ ] {r pl ,r p2 ,r p3 ) ^ 

v pi, P 2, P 3 w p (r pl 2)w p (r p 23)+w p (r p 23)w p (r p3 i)+Wp(r p31 )w p (r pl 2)' 



-8- 



So the relation between galaxy and matter reduced 3PCFs in equation (16) still holds, 

1 b 2 

Qp,9 = ~foQp,m + ' (19) 

if the reduced projected 3PCF is properly defined. This kind of definition is also adopted by Jing 
& Borner (1998, 2004), who point out that if Q is a constant (which in general it is not) then Q p 
is also a constant and equal to Q. Jing & Borner (1998, 2004) also measure the reduced projected 
3PCF as a function of triangle shape for galaxies in the LCRS and the 2dFGRS. 

Based on the above results, it is straightforward to infer bias from projected quantities. By 
projecting the measured redshift-space £ and £' 3 l, we form Q Pt9 for projected triangles that have 
two sides (r p i3 and r P 2z) fixed but differ in the angle $ between them. The dark matter Q P)m can be 
computed using the method introduced in § 2 (for projecting £ and t^) or can be measured from 
cosmological A^-body simulations. The reduced 3PCFs Q p , g {r p i3,r P 23,&) and Q p ,m(?pi3> t>23, $) 
can be expanded in either Fourier or multipole series, whichever one prefers. Their expansion 
coefficients satisfy Q p , 9 ,„(r p i 3 , r p23 ) = Q P ,m,n(r p i 3 , r p23 )/b + 5o n b 2 /b 2 , where 5 0n = 1 for n = 0, and 
otherwise. The nonlinear bias factor only enters in the n = component (i.e., the DC component 
for Fourier expansion or the monopole for multipole expansion). The relation between coefficients 
of other components provide useful tests of bias models and perturbation theory (see § 4 of Szapudi 
2004 for a more detailed discussion, which also applies here). Alternatively, we can infer b and b 2 
from a two-parameter fit to Q p>g (r p i3, r p2 3, $), which is widely used (e.g., Fry 1994; Feldman et al. 
2001; Scoccimarro et al. 2001). 

In order to make meaningful comparisons with predictions in the weakly nonlinear regime, all 
quantities should be measured on large scales. If the two sides (r p \z and ^23) of the projected trian- 
gle are large (e.g., greater than 10/i _1 Mpc) and differ substantially (e.g., |r p i3 — r P 23| > 5/i -1 Mpc), 
then the condition is guaranteed because any deprojected triangle lies in the weakly nonlinear 
regime (e.g., r > 5/i _1 Mpc). For configurations with r p i 3 « r p2 3 (but still large), the strongly 
nonlinear regime affects the value of Q p at <& ~ 0, where r p \ 2 approaches the size of the largest 
halos. Even for this extreme case, we can fit Q Pt g by rejecting data points near = 0. Since the 
strong nonlinear effect tends to add high-frequency components or higher order multipoles to Q p , 
it is unlikely for it to have a large impact on bias parameters estimated from Fourier coefficients or 
multipole coefficients. The use of projected triangles with different values of r p i 3 and r p23 would 
lead to a consistency check of the inferred bias factors and increase the signal-to-noise ratio. 

As an illustration, we fit the projected 3PCFs measured by Jing h Borner (2004) for galaxies 
in the 2dFGRS. Jing & Borner (2004) have measurements for two samples, one bright sample 
(Mf, < —19.66) and one faint sample (—19.66 < Mj, < 18.5). They also have predictions of reduced 
projected 3PCFs of the matter for the concordance cosmological model, which are measured in an 
A-body simulation. In their Figures 17 and 18, they compare the observed Q Pt9 and the predicted 
Qp,m for different triangle shapes, which provide the information we need to do the two-parameter 
fit. A different parameterization is adopted by Jing & Borner (2004) to characterize a triangle, one 
parameter (r p , the length of the shortest side) for the size and two parameters (u and v) for the 
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shape, with the lengths of the three sides being r p , ur p , and (u + v)r p . In their figures, the 3PCFs 
are plotted as a function of v (five equal linear bins in the range < v < 1) for several combinations 
of r p and u. For the two-parameter fit, we choose the case with the largest value of r p they have, 
i.e., 3.25/i _1 Mpc, and u = 2.09. If the v dependence is translated into an angular dependence, the 
five data points only cover $ from ~ 80° to ~ 180°, with the widths of the five angular bins being 
12°, 13°, 15°, 19°, and 44°. We see that the last bin smears the angular dependence a lot. So when 
measuring the 3PCFs from galaxy clustering data, adopting the (r p i3, r p 23, parameterization 
and dividing <I> into narrow bins are probably more suitable to probe the angular dependence than 
the (r p ,u, v) parameterization. Since the errorbars on the measurements are large and the scales 
are not truly in the weakly nonlinear regime, we cannot obtain robust constraints on galaxy bias 
parameters. Nevertheless, application of the formalism gives bias factors for galaxies in the bright 
(faint) sample that are consistent with b = 1.8 and 62 = (b = 1.1 and 62 = 0), somewhat higher b 
than expected but not absurdly so. 

If the projection is not infinite, the kind of relation shown in equation (16) also holds as long 
as we use the finite projected correlation functions in the definition of Q p (equation [18]). Strictly 
speaking, in this case, the integration of the product of 2PCFs in the right side of equation (17) 
cannot be written as the product of w p s except for ^(^13)^(^23)- However, the product of w p s 
should remain as a good approximation to the results. If one is not satisfied with the approximation, 
although it is good, one can always compute exact values of these integrations for Q Pt m and form 
Q P: g from the observation in the same way. Based on equation (17), we still have the relation in 
equation (19). 

The result can also be generalized to a projected field, where projected correlations include 
the effect of the selection function (see Fry & Thomas 1999). Based on equation (17), it is easy 
to show that equation (19) holds for Q p defined in terms of selection- function- weighted projected 
2PCFs and 3PCFs. Fry & Thomas (1999) perform a systematic study of projection effects on the 
reduced 3PCFs. They show that projections that are not deep enough would change the shape of 
the reduced 3PCF and thus bias the estimation of galaxy bias factors. However, this is based on the 
comparison with the three-dimensional reduced 3PCF of the matter. Our point here is that once 
the reduced 3PCF of the matter is calculated by taking account of the selection function, galaxy 
bias factors can be correctly inferred by comparing it with the observed reduced projected 3PCFs 
of galaxies. That is, we should always form Q Pi9 and Q p ^ m hi the same way. The depth of the survey 
should be much larger than the extent of structures caused by galaxy peculiar velocities so that 
we can calculate 3PCFs of the matter in real space. Buchalter, Kamionkowski, h Jaffe (2000) also 
investigate the reduced angular 3PCF. They concentrate on the explicit dependence of the galaxy 
3PCF on cosmological parameters, the selection function, and bias, not on the relationship between 
the reduced 3PCF of galaxies and that of the matter as done here. Under certain assumptions, 
their results can be put into the form discussed here. Focusing also on the explicit expression of 
the galaxy bispectrum, Verde, Heavens, & Matarrese (2000) present a theoretical analysis of the 
projected bispectrum in spherical harmonics and discuss its application (note that they find that 
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it is not encouraging to constrain bias parameters using projected galaxy catalogs). 

4. Summary and Discussion 

I investigate the projected 3PCF and find that it can be put in a simple form if Fourier 
expansion is introduced. Each Fourier component of the projected 3PCF is just a transform of 
the corresponding Fourier component of the bispectrum (equation [13]). In the weakly nonlinear 
regime, only the first few Fourier components are significant. The result is parallel to the multipole 
expansion proposed by Szapudi (2004) for the three-dimensional three-point statistics. Fourier 
expansion reduces the amount of work needed to compute the projected 3PCF and provides a 
convenient way to characterize its dependence on the triangle configuration. Projected 2PCFs 
and 3PCFs formed from redshift-space correlation functions are little affected by redshift-space 
distortions, which makes it promising to use them to measure galaxy bias factors. In the weakly 
nonlinear regime, I find that the relation between the reduced three-dimensional 3PCFs of galaxies 
and matter also holds for the reduced projected 3PCFs, if these are defined from the projected 
3PCF and 2PCF in the same way. The linear bias factor b and the first-order nonlinear bias factor 
&2 thus can be inferred from the reduced projected 3PCFs of galaxies (Q Pt g) formed from redshift- 
space measurements, through the dependence of Q P) g on the angle between two sides of projected 
triangles. From the point of view of Fourier expansion, the nonlinear bias factor only affects the 
DC component of Q p , and the linear bias factor is given by the ratio of coefficients of other Fourier 
components of Q p>g and Q p>m . 

The method of measuring galaxy bias factors from projected 3PCFs can be directly applied 
to the data from contemporary galaxy redshift surveys (e.g., 2dFGRS and SDSS). The projected 
3PCFs of galaxies should be measured on large scales, i.e., in the weakly nonlinear regime, where 
the bias relation (equation [19]) predicted by the local bias model applies. The projected 3PCFs 
of matter can be calculated using the technique introduced in § 2, or they can be measured from 
cosmological iV-body simulations. The inferred bias factors of galaxies would help to constrain 
cosmological parameters, such as the amplitude of the matter fluctuation power spectrum. 

Theoretically, the framework of the halo occupation distribution is often adopted to model 
galaxy clustering by linking galaxies to dark matter halos (for the three-dimensional 3PCF modeling 
within this framework, see, e.g., Takada & Jain 2003; Wang et al. 2004). With the help of this 
framework, even in the strongly non-linear regime, galaxy clustering would bring in additional 
constraining powers on cosmological parameters as well as on the relation between galaxies and 
matter. In any regime, projected 3PCFs, little affected by redshift-space distortions, are useful 
statistics. On very small scales where the three galaxies of each triplet come from the same dark 
matter halo, 3PCFs probe the shape of halos and the relative distribution of galaxies and matter 
inside halos. The projected 3PCF in this regime can be obtained by projecting the real-space 
3PCF, which is easy to calculate (Takada h, Jain 2003). On intermediate or large scales where the 
galaxies of each triplet are from two or three different halos, the calculation of the 3PCF is simpler 
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in Fourier space than in real space. Fourier expansion would be a useful technique to calculate the 
projected correlation functions on these scales. 

I am grateful to David Weinberg for his encouragement and advice on this work and for his 
detailed and valuable comments that greatly improved the paper. I thank Istvan Szapudi and Jun 
Pan for useful suggestions and comments. This work was supported by a Presidential Fellowship 
from the Graduate School of the Ohio State University and by NSF grant AST 04-07125. 




1.5 - 



-e- 

1— I 

m 



0.5 







5 
n 



10 




Fig. 1. — Weakly nonlinear matter bispectrum and Fourier expansion coefficients for k\ = 2A?2 = 
0.05/i _1 Mpc. The left panel shows the Fourier coefficients B n , which become virtually zero for 
n > 5. The right panel shows the bispectrum B (solid curve) as a function of the angle <\> between 
ki and k2 and bispectra reconstructed from the first few Fourier components. The bispectra from 
the Fourier series up to n = 3 (dashed curve), 5 (dotted curve), and 10 (dot-dashed curve, indis- 
tinguishable from the solid curve) have maximum fractional errors of ~7%, ~1.5%, and ~0.04%, 
respectively. A cold dark matter linear power spectrum is assumed (see the text). Both B and B n 
are normalized by dividing 3P(k±)P(k2) for display. 
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